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ABSTRACT 

Swift J 1626.6-5 156 is a Be/X-ray binary that was in outburst from December 2005 until November 2008. We 
have examined RXTEfPCA and HEXTE spectra of three long observations of this source taken early in its 
outburst, when the PCA 2-20 keV count rate was >70 counts s _1 PCU -1 , as well as several combined observa- 
tions from different stages of the outburst. The spectra are best fit with an absorbed cutoff power law with a 
~ 6.4 keV iron emission line and a Gaussian optical depth absorption line at -lOkeV We present strong evi- 
dence that this absorption-like feature is a cyclotron resonance scattering feature, making Swift J1626. 6-5156 
a new candidate cyclotron line source. The redshifted energy of ~10keV implies a magnetic field strength of 
~ 8.6(1 -Hz) x 10 G in the region of the accretion column close to the magnetic poles where the cyclotron line 
is produced. Analysis of phase averaged spectra spanning the duration of the outburst suggests a possible posi- 
tive correlation between the fundamental cyclotron energy and source luminosity. Phase resolved spectroscopy 
from a long observation reveals a variable cyclotron line energy, with phase dependence similar to a variety of 
other pulsars, as well as the first harmonic of the fundamental cyclotron line. 
Subject headings: stars: individual (Swift J1626. 6-5156) — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

The high mass X-ray binary (HMXB) Swift J1626.6-5156 
was discovered in outburst by th e Swift Burst Alert T elescope 
(BAT) on 18 December 2005 (iKrimm et al.l 120051) . It was 
soon recognized to be a n X-ray pulsar with a -15 s spin pe- 
riod dPalmer et alJl2005h and its compa nion was classified as 
a BOVe star at a dista nce of ~10kpc ( Negueruela & Marco 
l2006tlReig et al.ll2011l) . Shortly after the discovery, the Rossi 
X-ray Timing Explorer (RXTE) began monitoring the source 
and continued to do so for nearly five years. The catego- 
rization of Swift J1626. 6-5156 as a Be/X-ray (BeX) binary 
places it in a group that comprises most of the HMXB popula- 
tion, so we will describe it in terms associated with that group 
of objects. Its 2-20 keV PCA light curve is shown in FigureQ] 
spanning the beginning of RXTE's monitoring in 2006 Jan- 
uary through 2009 January, when the source had returned to 
quiescence. 

Swift J1626. 6-5156 was discovered during a Type II out- 
burst, which is thought to occur when the Be star's cir- 
cumstellar disk expands and temporarily engulfs the neutron 
star, leading to enhanced accretion and therefore a large in- 
crease in X-ray emission (Coe 2000, and references therein). 
iReig et all (120081) describe the exponential decay of the out- 
burst and the flaring behavior that occurred in the weeks 
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Figure 1. PCA light curve of Swift J1626.6-5156, taken 
from HEASAR C's RXTE mi ssion-long data archive 
I ftp://legacy.gsfc. nasa.gov/xte/data/archive/MissionLongDatai. It shows the 
background subtracted 2-20 keV count rate as recorded in standard2 mode 
(16 s time resolution), averaged over each observation and normalized to a 
single PCU. The three vertical lines indicate when the longest observations 
were taken in January 2006. LOl is the blue dashed line, L02 the green 
dot-dashed line, and L03 the red dotted line. The shaded regions represent 
the Data Sections, with DS1 in yellow, DS2 in orange, and DS3 in teal. 



immediately following the discovery. Additional flaring oc- 
curred following JD 2453800 (2006 March 5), after which the 
source transitioned into an epoch of long-term, quasi-periodic 
X-ray oscillations ( Reig et aT]j2008 ). A Lomb-Scargle peri- 
odogram ( Lomb ll976tlScarglel 1982) of the light curve first re- 
vealed an oscillation period from -JD 2454000-JD 2454350 
(2006 September-2007 September) of -47 days and from 
~JD 2454350-2454790 (2007 S eptember-2008 Nov ember) 
of -3/2 of t his, or -72.5 days dDeCesar et al.l |2009). Re- 



cent work bv lBavkal et al. 20101 found the long-term variabil- 
ity timescale to range from -45 to 95 days. Additionally, 
from pulse timing analysis, they find the binary orbit to be 
nearly circular (eccentricity 0.08) with a period of 132.9 days, 
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Table 1 

Observation Log 



Name 


Observation ID 


Date 


JD - 2450000 


PC A ? eX p 


HEXTE r CX p [ks] 


LOl 


91081-13-01 


2006 Jan 1 


3736.05 


26.4 


8.4 


L02 


91081-13-02 


2006 Jan 6 


3741.15 


23.6 


7.4 


DS1 


91082-01-01 -91082-01-27 


2006 Jan 12-30 


3747.73-3765.10 


23.3 


7.1 


L03 


91081-13-04 


2006 Jan 31 


3766.59 


27.6 


7.7 


DS2 


91082-01-32-91082-01-93 


2006 Feb 1-Mar 13 


3767.64-3807.52 


51.2 




DS3 


92412-01-27 - 93402-01-48 


2006 Jun 26-2008 Oct 10 


3912.21^1750.82 


78.8 





meaning the oscillations occur on timescales of ~l/2 and 
~2/3 of the orbital period. The oscillations therefore can- 
not be explained by bursts of accretion during periastron 
passages. On JD 2454791 (2008 November 11) the oscilla- 
tory behavior ceased and the 2-20 keV PC A count rate from 
Swift J1626.6-5156 dropped to ~1 counts 1 PCU" 1 . 

Many BeX and other HMXB systems harbor neutron 
stars with very strong magnetic fields, typical strengths be- 
ing ~ 10 n ~ 13 G, and some of these systems exhibit cy- 
clotron resonance scattering features (CRSFs or cyclotron 
lines). Observing cyclotron lines in the spectra of these 
systems is the only direct way to measure the magnetic 
field strengths of neutron stars. At the time of writ- 
ing there are sixt een neutron star X-ray binaries with 
confirmed CRSFs dHeindl et alj I2004t iCaballero & Wilms! 
120 12t iPottschmidt et al. 20121 and references therein), eight 
of which are transients (Swift J1626. 6-5156 being the 
ninth). Of the transients four systems are BeX binaries 
(SwiftJ1626. 6-5156 being the fifth) and four are OeX bina- 
ries. The persistent sources are HMXBs with an O- or B- 
star accompanying the neutron star (with exception of Her X- 
1 and 4U 1626-67, where the companion is of lower mass). 
The fundamental cyclotron line energy is related to the mag- 
netic field strength of the neutron star £ cyc ~ 11.6 Bi 2(1 + 
z)" 1 keV where B12 = fi/10 12 G and z is the surface gravita- 
tional redshift. Assuming a reasonable value of z — 0.3, corre- 
sponding to a typical neutron star mass of 1.4M , allows the 
direct calculation of B in the region local to the magnetic pole 
where the cyclotron lines are produced. 

In this Paper we discuss the discovery of a CRSF in the 
spectrum of Swift J1626. 6-5156. In |_2]we describe our data 
reduction and analysis techniques. $3] presents results from 
modeling the broad band X-ray spectra, including CRSFs 
and other spectral properties of the source. We summarize 
our findings and compare Swift J1626. 6-5156 with other cy- 
clotron line sources in ^4] and conclude in ^5] 

2. OBSERVATIONS AND DATA ANALYSIS 

For this analysis, we considered data f rom the Proportional 
Counter Array (PCA: |jahoda et al.l2006|) and the High Energy 
X-ray Timing Experiment (HEXTE: iRothschild et all 1 19981) 
onboard RXTE that were taken between December 2005 and 
April 2009. From HEXTE, we used spectra from only the B 
cluster, since the A cluster had stopped rocking at the time 
of these observation^ i.e., no background measurements for 
cluster A were available anymore. We extracted spectra in 
standard2 mode from the top layer of each PCA Propor- 
tional Counter Unit (PCU) separately. We then combined 
all spectra from a given observation into one, following the 

9 See http : //heasarc . gs£c . nasa . gov/docs/xte/whatsnew/big . html 
for more detailed information. 



recipe in the RXTE Cookbook and excluding data taken 
within lOmin of the South Atlantic Anomaly. For a more 
detailed de scription of the dat a reduction and screening pro- 
cedure see Wil ms etaTl ( 120061) . We considered the PCA en- 
ergy range 4-22 ke V and the HEXTE range 1 8-60 ke V when 
modeling the phase averaged spectra. 

Three long observations with PCA exposure time >20ks 
are available in the RXTE archive, along with hundreds of 
shorter monitoring observations. We combined the individual 
pointings of each observation as described in the RXTE Cook- 
book. The three "Long Observations", which we refer to as 
LO 1 , L02, and L03, were each analyzed separately, while the 
data between and beyond these observations were combined 
into "Data Sections", referred to as DS1, DS2, and DS3. Data 
prior to LOl were excluded in order to avoid spectral contam- 
ination from flar ing events (for a description of the flares, see 
iReig et all 120081) . For DS3 we combined only observations 
from the oscillating phase that had an average PCA 2-20 keV 
count rate of 10 counts s 1 PCU 1 or higher. The observation 
identification numbers and other information are given in Ta- 
ble Q] The locations in time of LOl-3 are shown by the col- 
ored lines and of DS1-3 by the shaded regions in Figure Q] 
The HEXTE data have low count rates, so we rebinned the 
counts at high energies using the FT00L qrppha^l. For LOl 
and L02 we rebinned by three energy channels above 50 keV 
and for L03 by four channels above 48 keV. We did not re- 
bin any counts for DSL The HEXTE data of DS2 and DS3 
were so noisy we decided not to use them as th ey might skew 
our spectral fits. We us ed xspecl2 models ( Arnaudl ll996t 
iDorman & Arnaudll2001l) for all spectral fitting 

3. SPECTRAL ANALYSIS 

3.1. Spectral Model 

We modeled the combined PCA and HEXTE (where 
present) phase averaged spectra of each LO and DS with 
a power law modified by a smooth exponential rollover 
(cutoffpl in xspec) including photoelectric absorption 
(phabs in xspec) and a Gaussian 6.4 keV Fe Ka emission 
line (gauss in xspec): 

/cont + Fe(£) = e^^ (a^/TV^ 

+ Ape e -j[(£-£Fc) 2 /<]| 
<x Fe V2tt J 

10 T he RXTE Cookbook can be found at 

http: //heasarc .gsfc. nasa. gov/docs/xte/recipes/cook_book. html 

11 See http://heasarc.nasa.gov/ftools/ftools_menu.html for 
more detailed information. 

12 The manual can be found at 

http : //heasarc . nasa . gov/xanadu/xspec/manual/manual . html 
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Figure 2. The ~4-60 keV spectrum and best-fit model for observation L03 
are shown in the main panel. PCA data are shown in red, HEXTE data in 
blue, a) The Ax 2 residuals left from fitting the spectrum with a cutoff power 
law. The spectrum peaks near 6.4 keV and dips near lOkeV. b) Residuals 
from adding a Gaussian emission line at ~ 6.4 keV to the cutoff power law, 
completing our continuum plus Fe line fit from Equation [T] The depression 
at 10 keV remains, c) Including a Gaussian optical depth absorption line near 
lOkeV removes the absorption-like feature and greatly improves the model, 
giving our best fit (Equation|2). There may be a second feature near 18keV. 
d) A second Gaussian optical depth absorption feature included near 18keV 
marginally improves the fit. The remai ning re sidual near ~ 36keV is most 
likely caused by an instrumental effect ( 33.2.71 . 



where Nn and are the absorption model components and 
are respectively the hydrogen column density per H atom for 
material of cosmic abundance (atoms cm -2 ) and the bound- 
free photoelectric a bsorption cross sec t ion (c m 2 ). We used 
the abund ances of [Anders & Grevessd d 19891) and th e cross 
sections of Balucihska-Church & McCammon (1992). 

The first term in Eq. Q]is the cutoff power law, where A cut 
is the power law normalization (photons keV -1 crrT 2 s _1 at 
1 keV), r is the power law photon index, and Etbid is the e- 
folding energy of the exponential rollover. The second term 
is the Gaussian emission line, in which is its normal- 
ization (photons crrT 2 s ), E Fe its energy, and cr Fe its width. 
This model was applied concurrently to the PCA and HEXTE 
data, differing only by a constant flux cross-calibration factor 
where the constant was fixed at 1 for the PCA and fit to ~0.8 
for HEXTE. For each spectral fit we also fit the background 
strength adjusting it on a level of ~0.1% using the command 
recornormin xspec. 

Residuals for a continuum fit using the model from Eq. Q] 
without and with the iron emission line for L03 are shown 
in Figs [2h and [2t>, respectively. A depression in the residu- 
als at ~10keV is apparent. In order to account for this fea- 
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Figure 3. Removing the Gaussian optical depth absorption component in 
our spectral model reveals the profile of the fundamental cyclotron line at 
~10keV in LOl (blue), L02 (green), and L03 (red). A decrease of the cen- 
tral energy between L01-L02 (close in time) and L03 is apparent, as well 
as an increase in line depth. 

ture, we modified the continuum plus Fe line model with a 
multiplicative absorption-like cyclotron line with a Gaussian 
optical depth profile (gabs in xspec): 



1(E) = 7 cont+Fe (£')exp 



l eye 



eye 



V2^ 



exp 



(E 



F ~) 2 

^cyc/ 



2a 2 

eye 



(2) 



Here r cyc , cr cyc , and E cyc are respectively the optical depth, 
line width, and line energy of the cyclotron line. The resid- 
uals of a fit to L03 using the full model 1(E) of Eq. [2] i.e., 
including the cyclotron line, are displayed in Figure [2]:. 

In the following sections we describe the modeling and 
interpretation of this line as a CRSF including several 
tests that led to the adoption of the model described 
above as the best-fit model for the phase averaged spectra 
(constarttxphabsx(cutof£pl + gauss)xgabs in xspec). 

3.2. Pulse phase averaged spectra 
3.2.1. Best-fit results, 10 keV CRSF detection 

We used L03, the longest observation of 
Swift J1626. 6-5156, to illustrate the results obtained 
with the best-fit model described in the previous section. 
This observation was taken on 2006 January 31, during the 
outburst decay (Figure [T]). The first PCA pointing is ~21 ks 
long and is followed immediately by a shorter pointing of 
~6ks. We combined them to form a single spectrum with 
PCA integration time ~27ks. The corresponding HEXTE 
observation has ~7.6ks of live time. The best-fit parameters 
for the L03 spectrum using the model of Eq. [2] along with 
the iron line equivalent width and other characteristics, can 
be found in the fourth row of Table [2] The phase averaged 
spectrum of L03 and best-fit model are shown in the main 
panel of Figure |2] Inspecting the residuals obtained after 
fitting the continuum model and iron line (EquationQ]) shown 
in Figure [2J3 more closely reveals the absorption-like feature 
near lOkeV and/or an emission-like feature near 14keV. 
Here we treated the feature as absorption-l ike only (Eq ua- 
tion |2} and obtained a very good fit. See {33. 2. 21kT2.5l for 
testing possible alternative descriptions. Fitting the feature 
with a Gaussian optical depth absorption line flattened the 
residuals in Figure |2j;, yielding x^. = 0.94 for 56 degrees 
of freedom. The common A^- 2 test available in xspec 
could not be applied in this case, as it is not applicab le to 
multiplicative model components (Orlandini et al.ll2012l) . We 
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Table 2 

Spectral parameters for phase averaged spectra 



Name TVh 


r 


^cut 


-Efold 








E C y C 


cr cyc 


1"cyc 


F x 2 /N iof 


[10 22 ]* 




[10 _1 ] c 


[keV] 


[keV] 


[keV] 


[10~Y 


[keV] 


[keV] 




[itr 10 cgs] 



LOl 


4 8 +0 - 6 
4 '°-0.6 


i i i +0.04 
1 ' 11 -0.05 


1 S+0-2 


10 46 +0 33 


6 33+0- 03 
°- JJ -0.04 


49+0- 05 

"■^-0.05 


T. 4+0- 5 
J,;, -0.4 


10 23 +ao '' 


99+0- 17 
"•"-0.17 






81.9/68 


L02 


H - J -0.7 


1 20 +0 - 06 


2 5 +03 

^-0.2 


11 50 +052 

11J "-0.39 


6 30 +0 - 04 

"■ J "-0.04 


54+0- 07 


3 +0 - 4 


10 01+ 09 
iu.ui _ 00g 


76 +015 

"'"-0.16 


0-082^ 


16.80!°;° 3 


78.9/68 


DS1 


4 8 +0 - 8 
-0.7 


i o 7 +0.07 
'■ -0.06 


9 4+0.3 
-0.3 


12 72 +0 - 79 
■ -0.37 


6 34 +0 - 05 

"-' -0.04 


45+ - 06 
-0.08 


1 7 +0 - 3 
'•'-0.3 


9 74+0.07 
y -0.07 


96+ 013 

"•^"-0.12 


15 +002 

"■ -0.02 


11.80$" 


94.1/76 


L03 


3 4 +0 - 7 
■'-0.8 


i qo+0.09 
1- -0.09 


1 7 +0 - 2 
'■'-0.2 


10 98+ 97 
iu.?o _ 09 , 


ft ^o+0.04 


19+ 010 




Q 71+0.07 
-0.07 


i 97+O.I6 
'■ -0.15 


A -3O+0.04 

"• -0.06 


o ft9 +0.04 
8.62 _ Q19 


52.8/56 


DS2 


1 4+0.8 
1 - -1.2 


1 61+ 007 

'■ ul -0.08 


1 o +(U 

i.u j 


19 07+1.83 

-2.47 


6 36 +0 - 03 

"■'"-0.06 


041 +0.05 
"■ -0.06 




10 00+ 08 


64+ 02 
" " -0.02 


25 +004 


4 44+0.03 
-0.15 


59.4/33 


DS3" 




i o 7 +0.07 
'•-"-0.07 


33+0- 02 

"■ J -0.02 


13 20 +1 67 


fi 4S+0.02 
6.48_ 002 


Q4I+0.05 


0.40^.04 


Q 44+0.07 
y - 44 -0.06 


0.64+°;'' 


25 +004 

"■ ZJ -0.04 


2 06 +0 003 

z -""-0.04 


52.4/32 



Note. — The best-fit spectral parameters from each observation or group of observations. is the column density, T the power law photon index, A cul the power law normalization, 
£V u i,] the folding energy of the cutoff power law, Ep c , o"f c , and Ap e the iron emission line energy, width, and normalization, Emc , <T cyc , and r cyc the energy, width, and optical depth of 
the Gaussian optical depth absorption line profile used to fit the cyclotron line, F the 3-20 keV source flux in units of 10~ 10 erg cm" 2 s" 1 , and^ 2 /A ? ( ] i the reduced^ 2 (where iVdof is 
the number of degrees of freedom), of each spectral fit. The changing A^r results from differences in rebinning of the HEXTE data (see for the L01-L03 columns, and from our 
choice to not use HEXTE data for DS2-3. The error bars quoted are at the 90% level. The 90% flux errors were calculated in xspec using 1000 draws from a Gaussian distribution. 
<7 DS3 requires an additional Gaussian emission component with energy 15.18+jj ]^ keV, tr fixed at 0.30 keV, and normalization 7.3+]^ X 10~ 5 photons cm" 2 s" 1 . * Values of are in 
units of 10 22 atoms cm" 2 . c Units are 10" 1 photons keV" 1 cm" 2 s" 1 at 1 keV. ''Values are given in units of 10" 3 photons cm" 2 s" 1 . 



instead used the more appropriate ratio of variances F-test 
dPress et alj|2007l) to calculate the probability that including 
this absorption-like featu re in the model i mprov es the fit by 
chance. As discussed by Orlandi ni et al.l (120121) . this F-test 
is not available in xspec; throughout this paper, we use 
the IDL routine mp£test(3 to calculate the probability of 
chance improvement (PCI) from one model to another. If one 
model is not a significant improvement over another, then the 
PCI will be large. The PCI of the model that includes the 
Gaussian-depth absorption feature (Figure [2};) compared to 
that without (Figure |2b) is 4.6 x 10~ 9 . The improvement to 
the fit due to use of the absorption-like model component is 
therefore very statistically significant. The best-fit parameters 
of this model are given in Table [2] 

The CRSF is detected in the earlier two long outburst ob- 
servations, LOl and L02, as well. Figure [3] shows the line 
profiles for L01-L03. We also detected a strong iron line in 
each of the spectra. See Table |2]for the best-fit results. One 
question we may ask is whether or not the cyclotron line per- 
sists throughout and beyond the primary outburst decay. To 
answer this we combined sections of data into single spectra 
to search for the line. The data selection is described in §2 
and summarized in Table Q] We found that both the iron line 
and the CRSF were present throughout the active phase of 
Swift J1626. 6-5156, as the features are detected in each of 
the Data Sections DS1-DS3 as well as in the Long Observa- 
tions. Interestingly, we found that the spectrum from the oscil- 
latory stage of the light curve (DS3) is best fit by the addition 
of a Gaussian emission line at ~ 15 keV, a component that was 
not necessar y in an y of the spectra taken during the outburst 
decay (see a3.2.5t bu t that is commonly used to fit HMXB 
spectra (Cobur n et al.l 120021) . All fit parameters are given in 
Table [2] The values of E cyc are suggestive of an evolution of 
the cycl otron line energy with luminosity, which is discussed 
in 34.31 Note that the parameters might be affected by ad- 
ditional uncertainties due to the range of fluxes covered for 
DS1-DS3 and missing HEXTE data for DS2-DS3, though. 

"Bumps" or "wiggles" in the spectra of accreting X-ray 
pulsars, especially near lOk eV, are common and have been 
discussed at some length by Cob urn et al.l (|2002). They find 
residuals in the data/model ratio at the ~0.8% level in all 
the accreting systems they examine. These residuals could 

13 http://www.physics.wisc . edu/$\sim$craigm/idl/down/mp£test 



easily be interpreted as cyclotron features. There have been 
marginal detections of absorption-like features, possibly cy- 
clotron lines, near lOkeV in other sources. One example is 
the rec ent work on XMMU J054 1 34.7-682550 by |InametiD 
(2009). These authors see what may be an absorption-like 
feature at lOkeV, but the count rate is too low to claim a sig- 
nificant detection. One must thus take good care to ensure 
that the feature we see is truly a cyclotron line and not due 
to calibration uncertainties, additional emission (e.g., from 
the Galactic ridge), or incorrect modeling of the continuum. 
We note that the residuals we see at ~10keV are on the 3- 
10% level (Figure[3]l, i.e., comparatively strong. Nevert heless 
we add ress t hese po t ential sources of uncertainties in 33.2.21 
31231 and 3T2~4l - 3T231 respectively. We also note that 
Reig e t al.l (1201 11) state that they do not detect a cyclotron line 
in the RXTE outburst data. However, these authors include 
a 9keV edge in their model apparently masking the CRSF. 
This model choice misses the typical CRSF properties of the 
observed residual that we present in this paper (e.g., its lumi- 
nosity or pulse phase dependence). 

3.2.2. Crab spectrum comparison 

One test we can perform in order to validate the existence 
of the cyclotron line in the phase averaged spectrum is to com- 
pare the spectrum of Swift J1626. 6-5156 with that of a non- 
cyclotron line source like the Crab pulsar. To be as consistent 
as possible in our comparison, we used Swift J1626. 6-5156 
and Crab data on ly from PCU 2, wh ich is the best calibrated 
of the five PCUs dJahoda et alJl2006T) and which is most often 
turned on during observations. The Crab data included both 
the pulsar and the nebula, and were taken between November 
2005 and March 2006, the same time period during which the 
Swift J1626. 6-5156 observations LOl-3 were taken. 

We first searched for an absorption-like feature near lOkeV 
in the Crab spectral residuals. The Crab has no CRSFs, so 
if such a feature were seen in the Crab spectrum, we could 
assume that the Swift J 1626.6-5 156 feature was instrumental 
in origin. We modeled the Crab spectrum with an absorbed 
power law, freez ing A^h at 4 x 10 2 1 cm" 2 , the value determined 
for the Crab by Weisskopf et al. (2004). The power law pa- 
rameters were left free, and we found that F = 2.095 with 
a normalization of 10.343 photons keV" 1 cm" 2 s" 1 at 1 keV 
gave the best fit. The data/model ratio is ~0.99 near lOkeV, 
P r ®hich is consistent with expected calibration uncertainties 
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Figure 4. Ratio of Swift J1626-515 6 (observation L03) to Crab counts ac- 
cording to the procedure described by Orlandini et al. ( 1998). The depression 
at ~10keV is present in this model- and calibration-independent representa- 
tion of the spectrum of Swift J1626-5156 as well. 



Table 3 

Comparison with alternative continuum models 



Continuum Model 




E C yc 


Y 2 


AW 


PCI 




[keV] 


[keV] 






cutoffpl 


ft ^o+0.04 
638 -0.04 


q 7 1 +0.07 
-0.07 


0.94 


56 




powerxhighecut 


ft oc+0.04 


q 74+O.O8 
'■'^-0.07 


0.96 


55 


0.47 


power xfdcut 


ft OO+0.04 


q 70+0.11 
y - -0.12 


1.12 


55 


0.26 


npex 


ft oq+0.06 


9 83+ - 28 


1.02 


54 


0.38 



Note. — Comparison between the cutoffpl continuum model with the 
highecut, fdcut, and npex models described in j l3.2.4l All spectral mod- 
eling was performed on the L03 dataset. Fits with the alternative continuum 
models yield values of £pe and Eqvc , defined as in Tableful that are consistent 
with the line energies found using the cutoffpl model. The probability of 
chance improvement (PCI) is described in j !3.2.4l and calculated with x^. and 
A^or of the cutoffpl and each of the alternative models. The PCI is large for 
each of the three alternative models, meaning that the different continuum mod- 
els are statistically equivalent. For ease of comparison, we repeat the selected 
results for the cutoffpl continuum from Table[2l 



( Jahod a et alj|2006h . No absorption-like feature is found at 
lOkeV in the Crab spectrum (and no rollover). We repeated 
the procedure with PCU 2 data from Swift J1626. 6-5156 dur- 
ing L03, using a cutoff power law and fixing A^h and F to the 
values in Table [2] In this case the ratio near lOkeV is ~0.91, 
a much larger residual than consistent with calibration errors. 

We next took a model and calibration independent approach 
by comparing the raw counts of Swift J1626. 6-5156 with 
those of the Crab in each spectral energy bin following the 
procedure of lOrlandini et al.l (119981) . We first subtracted the 
background counts from each raw spectrum, making sure to 
account for the difference in exposure time between the two 
sources. We then divided the Swift J1626. 6-5156 count rates 
by those of the Crab, resulting in the ratio shown in Fig- 
ure |U At <6keV the J1626/Crab ratio reflects the steeper 
slope of the Crab and at >14keV the missing Crab rollover. 
Between these energies an enhancement at 6.4 keV due to 
the iron line in Swift J1626. 6-5156 and a clear depression 
at ~10keV, rather than a smooth transition between the iron 
line and the power law tail, are apparent. The same kind of 



Crab ratio is shown, e.g., for the -55 keV CRSF of Ve la X-l 
in the upper panel of Figure 3 of lOrlandini et al.l d 1998b . This 
plot shows similar cyclotron line and rollover features as the 
J1626/Crab ratio (it differs at lower energies due to the overall 
different energy range). Independently of the response func- 
tion or spectral model the J1626/Crab ratio therefore supports 
that the absorption-like feature does not result from a calibra- 
tion error. 

3.2.3. Galactic ridge emission 

We also investigated the possibility that cyclotron features 
in the spectrum arise from the Galactic ridge emission by 
analyzing data from when the source was quiescent. While 
Swift J1626. 6-2126 remained fairly bright and highly vari- 
able for about 1000 days after the initial outburst, it then 
quickly faded to quiescence (Figure[TJ and stayed reliably low 
for the last two years of the RXTE mission. During this time 
the source continued to be monitored regularly by the PCA. 
We combined these late observations to obtain a spectrum 
with a total of 42 ks exposure, and compared this with our out- 
burst observations, in order to put limits on the contribution of 
the flux during the outburst observations from the underlying 
Galactic ridge emission. In the lOkeV region surrounding the 
fundamental cyclotron line, the ridge emission has a flux of 
less than 1% of the count rate of Swift 1626.6-5156 during 
our observations (L03). In the energ y ran ge around the first 
harmonic (~18keV, see 33.2.61 and 33. 3t no emission from 
the ridge is detected. Therefore Galactic ridge emission can- 
not explain the deviations that we measure. 

3.2.4. Alternative standard continuum and CRSF models 

In order to test the continuum model dependence of our 
results we in turn replaced the cutoff power law model by 
three other standard empirical pulsar continuum models and 
repeated the fit for L03. The first of these was the power law 
with a high energy cutoff (powerxhighecut in xspec where, 
in contrast to the smooth cutoff power law, the rollover in 
the s pectrum only starts at a certain energy, £ cut dWhite et al.l 
119831) . While this model is often used, it can create an arti- 
ficial line-li ke residual near E„ n due to the discontinuity at 
that energy (ICoburn et al.ll2002L and references therei n). Next 
we fit the spectrum with the Fermi-Dirac cutoff from Tanaka 
( 1986) (powerxfdcut in xspec). This model has a smooth, 
continuous rollover described by E cat and fit'dd, which are, 
however, not directly comparable to the values obtained with 
cutoffpl or powerxhighecut due to differences in the con- 
tinua. The final mod el we used is the Negative Positive Ex- 
ponential model from Mi haral d 19951) (npex in xspec), which 
consists of two power laws with a smooth exponential cut- 
off described by Zifoid- Each of these models was fit to L03 
modified by absorption and with the added 6.4 keV Gaussian 
emission line and cyclotron line. The fit values of the Gaus- 
sian emission line and the CRSF energy found for each model 
are given in Table [3] 

These alternative continuum models were found to pro- 
duce fits of similar quality when compared with each other 
and with the cutoffpl model. They also produced identi- 
cal results, within errors, for all comparable parameters, most 
notably those of the lOkeV CRSF. The Ay 2 test cannot be 
used t o compare different continuum models (Protass ov et al.l 
120021: lOrlandini etaT]|2012l) . To quantify whether or not the 
cutoffpl model, which has the lowest ^ 2 ed of the four contin- 
uum models tested, actually improves the continuum fit when 
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compa red to the alternative models, we calculated the PCI 
(as in ED of the cutof fpl model. lOrlandini etail d2012l) 
discuss the appropriateness of this F-test for comparison be- 
tw een different models , and the statistic is applied in this way 
bv llwakiri et al.l (120121) . Table[3]gives the reduced x^ ed , num- 
ber of degrees of freedom, and PCI (compared to cutof fpl) 
for each alternative model. We find that the PCI is between 
~25-50%, meaning that the cutoffpl improves the contin- 
uum fit marginally at best. Because none of the canonical 
models clearly gave a best fit to the continuum, nor adjusted 
the energies of the 6.4keV line or the lOkeV feature, and 
since cutoffpl has the smallest number of parameters, we 
defaulted to this model, our original choice, as a result of this 
test. The cutoffpl model is also well suited for comparisons 
with many earlier results for similar sour ces ( Cobu rn et al. 
2002 1 ; iMowlavi et alj|2006t iCaballero et alj|2008t iSuchv et al. 
20ll . 



4 -(a), 



Also, in order to further check that the lOkeV feature is 
robustly described as a typical CRSF, we compared the Gaus- 
sian optical depth absorptio n model with a Lorentzian profile 
model (cyclabs in xspec) (Mihar aTl990l) . Taking the known 
systematic shift into account we found similar centroid ener- 
gies for the absorption-like feature using either model. We 
chose to use gabs in our analysis because its parametriza- 
tion provides an easier to interpret and compare characteris- 
tic energy parameter (in contrast to gabs the energy parame- 
ter of cyclabs does not directly reflect the energy where the 
absorption-like feature in the spectrum is deepest). 

3.2.5. Modified continuum models 

Many accreting pulsar spectra are either not adequately 
or at least not uniquely described by the standard empirical 
continuum models above, especially in the 10-20keV range 
dCoburn et al.H2002t lOrlandinill2004h . Incorrect modeling of 
the exact shape of the rollover can create residuals in addition 
to CRSF features or, in the worst case, mimic CRSFs. There 
are many examples for the application of con tinuum models 
with s lightly modified cutoff characterization. Orlandini et al. 
d 19991) noted the likely presence of a two step change of 
the spectral slope for OAO 1657-415 (as well as a poten- 
tial cyclotron line residual near 36 keV), the first change oc- 
curring at 10-20keV and the other one at higher energies, 
however, they were able to model this behavior taking ad- 
vantage of the discontinuity and the two characteristic ener- 
gies of the highecut model mentioned above. Other au- 
thors used an additional polynomial compone nt in order to 
describe a more complex rollover shape, i.e., iBurderi et al.l 
(2000) for Cen X-3 (an e stablished cyclotron line source) 
and lKlochkov etal.1 d2008l) for EXO 2030+375 (noting a po- 
tential cyclotron line residual near 63keV). The correct con- 
tinuum and therefore the presence of a 25keV CRSF in 
Vela X-l (in additio n to the established ~55keV one) has 
long been debated dOrlandini et alJfl998l : iRreyken bohm et al. 
20021 lLa Barbera et alj 120031) . In the most recent analy- 
sis its presence was co nfirmed in a Suzaku observation 
dDoroshenko et al.l 1201 ll) . Another widely used modifica- 
tion is obtained by adding a broad Gaussian component, 
most often in emission, to a standard or already modi- 
fied model dCoburn et all 120021: iKlochkoy et al.l l2007l l2008t 
iFerrigno et al l 12009b ISuchv et all 120081 1201 ll) . Since our 
Swift J1626. 6-5156 observations show a positive residual at 
~14keV before modeling the cyclotron line (Figure [2J3) we 
applied this modified model to test whether our cyclotron 
line detection is robust against the choice of a more complex 
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Figure 5. a) Starting from the continuum plus Fe line fit (Equation [1} for 
observation L03 (residuals shown in Figure [2p), we further smoothed the 
continuum in this alternative fit by first including a broad Gaussian emission 
feature centered at 13.7 keV. While this flattens the residuals substantially, the 
absorption-like feature at 10 keV still remains, b) Final fit for this mixed emis- 
sion/absorption approach after also including the 9.7 keV cyclotron line. The 
remaini ng res idual near ~ 36 keV is most likely caused by an instrumental 
effect ( j}XT7l . 



rollover description. 

We therefore altered our fit procedure by first adding a 
broad emission line near 14 keV to the initial continuum plus 
Fe line fit for L03 (Figure [2J3, Equation [TJ. The residuals 
of this fit are shown in Figure [5^. The negative residual 
near lOkeV is still clearly visible, giving a^ ed of 2.46 for 
56 dof. Thus an absorption-like feature was still required 
for an adequate description. Multiplying by a Gaussian op- 
tical depth absorption feature and refitting we obtained E cyc = 



9.74+°;}*keV, cr cyc 



0.93!°;gkeV, 



• eye 



for the 

' +0.981 



[ 6 keV, cr em = 1.86^keV, 



3 photons cm 2 s 1 for the emission-like 



cyclotron line and E em - 13.8 +0 

"em — U - Z ' <J -0.15 

line. This fit with both lines yielded x z /N do f = 49.96/53 
(= 0.94) (Figure Et), as opposed to 52.76/56 (= 0.94) for 
the original best-fit (Figure [2}:, Table The fits are equally 
good (PCI = 50%). While the depth of the cyclotron line was 
smaller in the modified fit, the line was still detected at the 
same energy as before (Table|2]i. The parameters of the modi- 
fied fit characterizing the continuum shape were F = 1.47 +q°j 
and Efoid = ^.SS^HfjkeV, consistent within errors with the 
original fit (Table |2). The same Fe line parameters were ob- 
tained for both fit procedures. 

This implies that, while the underlying continuum might 
be slightly different from a standard cutoff power law, such a 
difference cannot be significantly detected in this observation. 
The absorption-like feature at lOkeV on the other hand is sig- 
nificant. It appears at 9.7 keV regardless of the presence or 
absence of a broad emission feature and its energy is stable. 
Based on these results and since the model without the broad 
emission feature is the simpler one, it is the one we chose as 
our best-fit model. The same overall behavior has been seen 
for all the datasets with the exception of DS3, for which a 
better fit was obtained including a broad emission feature, see 
notes of Table [2] for further details. This can be understood 
since the DS3 spectrum was obtained averaging over many 
datasets, which might contain some evolution of the contin- 
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3.2.6. 1 8 keV harmonic 

Returning to Figure[2j: we point out a small depression near 
18keV. In Figure |2}i we show the obtained residuals after fit- 
ting this feature with another Gaussian optical depth absorp- 
tion line with centroid energy 18.5^?keV, r ~ 0.1too7> anc ^ 
<x, which could not be constrained, fixed at ~0.02. This addi- 
tional component reduced^ 2 slightly but did not significantly 
improv e the fit, as its PC I ~ 40% from the ratio of variances 
F-test (Pr ess et al.l 120071) . We therefore did not include this 
feature in our best-fit model of the phase averaged spectrum, 
i.e., Table|2] While we would not normally consider this a de- 
tection, we also found that the addition of this 18 keV feature 
greatly improved the pulse phase resolved spectral fits ( 83. 3b 
during the pulse peak. It is marginally present in the phase 
averaged spectrum of L02 as well. It is thus likely the first 
harmonic of the fundamental CRSF, which i s sig nificantly de- 
tected in only a few profile phase bins (see 83.3b . 

3.2.7. 36 keV instrumental feature 

Even after fitting the harmonic feature there remained an 
even more obvious line-like residual near 36keV (Figures |2jl 
and [St). It is tempting to interpret this as a further harmonic. 
There might even be another weak residual at <30 keV. How- 
ever, the effective area of HEXTE shows a sharp drop at 
33 keV due to the K-edge of iodine with a smooth recovery 
until about 45 keV(see Figure 5 of iRothschild et al.l 119981) . 
So there is a good chance that (at least part of) any resid- 
uals in this energy rang e (is) are due to imperfect calibra- 
tion dHeindl et al.lll999ab . In addition the strongest internal 
HEXTE background line, which is due to K-lines from the 
tellurium daughters of various iodine decays, sits at 30 keV 
(IRothschild et al.lll998l) . We therefore decided to not further 
pursue the potential existence of harmonics beyond the first 
one at 18keV. If Swift J1626. 5-5156 shows other outbursts 
in the future the detection of these harmonics might be possi- 
ble with higher sensitivity with Suzaku, or better, NuSTAR or 
Astro-H. 

3.3. Pulse phase resolved spectra 

We again used L03, which shows the strongest cyclotron 
line, for a phase resolved spectral analysis. We extracted 
Good Xenon event data with a time resolution of 3.9 ms and 
corrected the event times to the Solar System barycenter. A 
period search on these data resulted in a pulse period of 
15.35828 ± 0. 00008 s, consisten t with the spin down evolution 
presented by IBavkal et al. 20101 We folded the 3.5-25 keV 
events on this period in order to obtain the pulse profile, see 
dashed lin e in th e upper panel of Figure [6] As presented by 
iReig et all d2008l) the pulse profile shows a single broad peak 
with a roughly constant maximum over four phase bins. The 
falling edge seems to be steeper than the rising edge, which 
may be an artifact from the period folding. We did not in- 
clude a binary cor rection since effects of the ~133d orbit 
(IBavkal et al. 20101) are negligible for our <30 ks observation. 

We divided the pulse profile into 16 phase bins for spec- 
tral analysis and used ikf asebin to extract spectra for each 
phas^3. Each spectrum was fit using Eq.|2]with an additional 
edge component at ~5 keV to account for the Xenon L-edge 
calibration feature (Jahoda et al. 2006). Still we found it nec- 
essary to remove one spectral bin at ~4.5 keV (spectral bin 12) 

14 See http://pulsar.sternwarte.uni-erlangen.de/wilms/research/analysis/ 
rxte/pulse.html for more detailed information. 




Pulse Phase 



Figure 6. Variability of the spectral parameters with phase for obser- 
vation L03. The top panel includes the 3.5-25 keV pulse profile of 
Swift J1626. 6-5156 (dashed line) as well as the best-fit power law index 
(crosses). The shaded region indicates the four phase bins surrounding the 
pulse peak, in which a second CRSF line is necessary to produce an accept- 
able fit. 

from each phase's spectrum to improve the overall fit. We as- 
sume that the systematic residual in this spectral bin is due to 
calibration uncertainties at lower energies and has no physical 
significance for the model. The column density and the iron 
emission energy were frozen at 3.4 x 10 22 crrT 2 and 6.4keV - 
the values from the phase averaged fit -, respectively. In order 
to better constrain the spectral index, the folding energy of the 
cutoff power law was also frozen at the phase averaged value 
of 10.98 keV. 

Using only the continuum plus Fe line model (Eq. [TJ on 
the phase resolved spectra led to unacceptable fits throughout 
the whole pulse, see column ^ 2 d () in Table [4] Including a 
CRSF in the form of a Gaussian optical depth absorption line 
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Table 4 

Spectral parameters for phase resolved spectra from L03 



Phase 




r 




^cut 


crpc 




^cyc.l 


C cyc.l 




£cyc,2 


C cyc,2 


' cyc,2 


*red,2 


^red.l 


*red,0 


PCI'' 




Bin 






[lO-'f 


[keV] 


[lcr 3 ] 6 


[keV] 


[keV] 




[keV] 


[keV] 




41 dof 


43 dof 


46 dof 






1 


1 


14+0.02 
J -0.01 


1 


09+0.03 

■ -0.03 


3 +(u 


1 02+ - 50 
1 - uz -0.34 


10 2 +03 
lu - -0.2 


1 1+0.4 
1 - J -0.5 


4 +0 - 2 
-0.1 










1.34 


2.76 


9.2 x 10" 


-3 


2 


1 


18+0.02 
' -0.01 


1 


99+0.03 
-0.03 


3 +0 - 2 


1 15+ - 75 


10 0+ 03 
iu.u_ 02 


1 6 +0 5 


5 +01 










0.98 


2.10 


6.5 x 10" 


-3 


3 


1 


49+O.OI 
-0.01 


1 


57+0.03 
-0.03 


3 +0 " 2 


1 19+O.69 
1J -0.38 


9 6 +0 - 2 
y -°-o.i 


1 4+O.6 


5 +0 2 










0.93 


2.41 


1.0 x 10" 


-3 


4 


1 


•jo +0.01 


1 


1 +0.04 
■° -0.03 


+0 ° 


1 10 +OM 
1-1 -0.20 


9 6 +al 




5 +0 ° 










0.68 


4.64 


1.8 x 10" 


-9 


5 


1 


•29+O.OI 
' -0.01 


1 


09+O.O3 
-0 03 


°« 


1 15+ - 26 


9 8 +ai 
-0.1 


°m 


3+ 00 
u - -0.0 










0.96 


6.57 


1.7 x 10" 


-9 


6 


1 


90+0.01 

' -0.01 


1 


Si+0.04 
-0.03 




9 7C+O.52 
/J -0.47 


10 o +0 - 2 


4 +0 2 
-0.1 


2 +0 
u - z -o.o 










1.69 


5.46 


8.6 x 10" 


-9 


7 


1 


90+0.00 

' 28 -0.01 


1 


gn+0.03 


1 +ai 


2 9 6+o.5i 

^■'"-0.46 


10 +ai 


3 +0 - 2 


2 +0 ° 










1.41 


4.35 


1.5 x 10" 


-4 


8 


1 


9^+0.01 
-0.02 


1 


04+O.O3 
■ S4 -0.04 


o +0 - 4 


62+ - 26 

u u -0.19 


7+0.1 
'■'-0.1 


1 +0 2 
1,u -0.1 


4 +o.i 
u - -0.6 


17 5 +0 ' 6 

1 '- J -0.4 


2.5 




1.24 


2.29 


5.38 


2.5 x 10" 


-2 


9 


1 


99+0.02 
• -0.01 


1 


os+0.04 
-0.04 


0+ 03 


74+ 49 
-0.20 


1+0.2 
"-0.1 


1 3 +a2 
*• -0.2 


5 +1 ° 
-0.1 


18.4±*} 


2.5 




1.05 


2.80 


5.42 


1.0 x 10" 


-3 


10 


1 


99+0.03 
• -0.3 


1 


oc+0.07 
■° -0.08 


4 +al 


9 90+1. 14 
2.28_ 039 


9 8 +0 - 4 
'■°-0.4 


2 j+0.7 






2.5 


1 3+ 03 

i - : '-0.3 


1.33 


2.81 


3.64 


8.8 x 10" 


-3 


11 


1 


9 n +0.03 
-0.03 


1 


70+O.O8 
■ /o -0.08 


4 +ai 
u -^-o.o 


4 19+I.I6 
-1.02 


as*! 


9 9+0.6 
■ -0.5 


0.8!gj 


19.133 


2.5 


1 5 +0 - 4 


1.17 


2.17 


2.62 


2.5 x 10" 


-2 


12 


1 


9O+0.01 
'-0.01 


1 


Q4+0.03 
' -0.03 


4 +al 
u - -0.0 


4 40+0.68 
^ w -0.71 


11.8*0.3 


1 9+0.4 

-0.3 


2 +01 

-0.0 










1.43 


2.58 


2.7 x 10" 


-2 


13 


1 


90+O.OI 
' iS -0.01 


1 


fi o+0.03 
■ bS -0.03 


4 +ai 
u - -0.1 


2 44+0.85 
^ -0.69 


11 5+ 03 

11 •-'-0.4 


1 7 +a7 
1 ■ ' -0.4 


4 +0 2 
"• -0.1 










0.89 


2.25 


1.3 x 10" 


-3 


14 


1 


26 +o.oi 
• -0.01 


1 


1^+0.03 
J -0.03 




1 13+ - 71 

1 - 1 -'-0.34 


10 9+ 03 
lu - -0.3 


1 S+0-4 
L8 -0.4 


6 +0 2 

u '°-0.1 










1.03 


3.66 


2.6 x 10" 


-5 


15 


1 


•31+0.01 

-0.01 


1 


1 r+0.02 
■ -0.03 


3 +a2 


1 ic+0.45 
1.18_ 027 


10 5 +a2 




4 +o.i 










1.02 


3.49 


4.2 x 10" 


-5 


16 


1 


14+0.02 
' J -0.01 


1 


Oq+0.03 
■ u -0.02 


3 +a2 

""'-0.1 


1 1 1 +0.46 
^ -0.34 


10 2 +a2 
lu - -0.2 


0.8^ 


3 +01 










0.99 


2.38 


2.2 x 10" 


-3 



Note. — A C[ll and are respectively the normalization constants of the cutoff power law and Gaussian iron emission line; E Lyc \ , cr™ c> i , and T cyc j the energy, width, and Gaussian 
optical depth of the first CRSF; £ C yc,2, 0-cyc,2, and r C y C .2 the same parameters for the second CRSF; ^ re d,o, .Yred.l, and;\f m i,2 the reduced^ 2 values for fits with no CRSF, one (the 
fundamental) CRSF, and both CRSFs. "Units are 10~' photons keV -1 cm" 2 at 1 keV. ''Units are 10~ 3 total photons cm -2 s~' in the line. N^, £pe. ana ^lold were frozen to their phase 
averaged values. c As stated in the text in j !3.3l the probability of chance improvement (PCI) is calculated between the models with no CRSF and one CRSF for phase bins 1-7 and 
12-16, and between one CRSF and two CRSFs for phase bins 8-1 1 (the latter are listed in bol d print ). Th e PCI values between models with no CRSF and one CRSF for bins 8, 9, 10, 
and 11 are respectively 2.8 X 1(T 3 , 1.5 X 1CT 2 , 0.20, and 0.27. The PCI is described briefly in ETTI and j f3~2~4l 



at * 10keV (Eq. 13 improved the fits significantly, but did 
not result in a sufficiently low x 2 during the maximum of the 
main peak, see column xLa 1 m Table U Including a second 
absorption line in these four phase bins improved the fits to 
acceptable values, see column x 2 ed 2 m Table [4] The results 
from these fits are shown in Figure [6] and Table |4] The PCI 
dPress et al.l 120071) in the last column was calculated for the 
best fit in each phase bin, i.e., between the models with no 
CRSF and one CRSF for phase bins 1-7 and 12-16, and be- 
tween no CRSF and two CRSFs for phase bins 8-11. The four 
phase bins requiring the second cyclotron line are highlighted 
by the shaded region in Figure[6] The width of the second line 
was frozen at 2.5 keV in order to better constrain the other pa- 
rameters. The centroid energy is roughly twice the energy 
of the fundamental line and is consistent with being its first 
harmonic. The phase resolved results for the harmonic likely 
explain the additional feature mention ed in the description of 
the phase averaged spectrum ( 33.2.61 Figure |2jl). It has such 
a low significance in the phase averaged spectrum because it 
is not detected at most rotation phases and its energy changes 
with phase. There are other CRSFs that are only detected in 
phase resolved spectra. For example, the second harmonic of 
the CRSF in Her X-l is detected only in the des cent of the 
main peak dDi Salvo et al.ll2004tlEn"oto et al.H2008l) . 

Figure [6] shows that the power law index hardens through- 
out the peak and softens in the minimum. As can be expected 
the power law normalization roughly follows the pulse profile, 
reaching ~0. 18 photons cm _2 keV _1 in the peak. The CRSF 
energy shows a smaller increase during the rising edge, then 
first a decline and then a drastic increase of more than 20% 
in the peak itself, followed by a slow decline over the falling 



edge and minimum of the pulse profile. The width and depth 
of the CRSF are variable as well with a dip in the rising edge 
and a possible increase through the peak and a rapid decline 
at the beginning of the falling edge. 

4. DISCUSSION 

4. 1 . Magnetic field estimated from E cyc 

We found evidence for cyclotron resonance scattering fea- 
tures at ~10 and 18keV in the spectrum of the Be/X-ray bi- 
nary Swift J 1626.6-5 156. While the harmonic CRSF was not 
seen at a significant level in the phase averaged spectrum, 
pulse phase resolved spectroscopy showed it is a necessary 
addition to the spectral fit between phases ~0.4-0.7, i.e., dur- 
ing the pulse peak. We can use the cyclotron line energy to 
estimate the strength of the magnetic field local to the line 
production region near the neutron star's polar regions. The 
fundamental and subsequent harmonic cyclotron lin e energies 
are re lated to the local magnetic field strength by (M eszarosl 
1992) 



2 V 1+2 "(B/BLni)sin 2 fl-1 1 



1+Z 



(3) 



«(11.6keV)(l +z)- 1 Bi2 



where m e is the electron rest mass, n = (1, 2, 3, . . .) is the in- 
teger harmonic number, B cl , t ~4xl0 13 Gis the critical field 
strength for resonant scattering, is the angle between the 
photon direction and the magnetic field vector, z is the gravi- 
tational redshift at the neutron star surface (z ~ 0.3 for a neu- 
tron star with mass 1.4 M and radius ~ 10 6 cm) and B\2 = 
B/(10 12 G). A fundamental energy of ~10keV therefore im- 
plies a magnetic field of ~ 8.6(1 + z)10 n G ~ 1.1 x 10 12 G. 
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4.2. Magnetic field estimated from accretion 

As a consistency check we estimated the global magnetic 
field strength, assuming that the neutron star has reached 
its equilibrium spin rate through interplay b etween an accre- 
tion di sk a nd the magnetic field. We ref er to iGhosh & Lambl 
(fl979l) and IShapiro & Teukolskvl (fl98l Chapter 15) for our 
calculation and use the pulse period of L03. While the 
pulse period evolved over the outburst, the change was small 
enough to not throw the spin o ut of equilibrium: we refer to 
Figure 3 of Ba vkal et al. 2 010, which shows a change of only 
0.2% in the pulsar's spin period over the whole activity phase. 
This small change in P sp i n had little effect on our estimate of 
the global magnetic field. We estimated B e i bai by equating 
the Alfven radius r& with the corotation radius at which the 
orbital velocity is equivalent to the surface velocity of the ro- 
tating star, thus giving w sp ; n * ^GM/r 3 A rads~' where 

r A ^3.2xlO 8 M- 7 2/7 4 7 (0 1/7 cm (4) 

Here M [7 = M/(10 I7 gs~') is the mass accretion rate and 
yU30 = ju/ ( 10 30 Gem 3 ) is the dipolar magnetic moment. P sp j n = 
15.35828 s is our measured spin period, so <y sp j n ~ 0.41 rad s 
and r\ ~ 1 x 10 9 cm. To estimate the average M over the 
source's lifetime, we must firs t estimate the lum inosity of 
the source outside of outburst. iReig et all (120111) estimated 
a source distance of d = 10.7 + 3.5 kpc. We assume that the 
current low luminosity state of the source is typical outside 
of outburst and from the source's spectrum during this state 
we derived a 2-60 keV flux of 2.9 x 10~ n erg cirT 2 s _1 , giv- 
ing a luminosity of ~ 1 .7 x 10 35 erg s _1 at 10 kpc for isotropic 
emission. The mass accretion rate for an efficiency rj w 0.1 
is then M - L/(i]c 2 ) ~ 2 x 10 15 gs~'. Using Equation @] we 
find ^30 ~ 1 and the global dipolar magnetic field strength 
-Bgiobai - 2yu/-Ri s = 2 x 10 12 G. As a rough estimate this is rea- 
sonably close to the field strength derived from the cyclotron 
line (see previous section). 

4.3. E cyc dependence on luminosity 

Something else to consider is the relationship between 
the cyclotron line energy and the source's X-ray luminos- 
ity. Her X-l and GX 304-1 are examples of systems 
containing an accreting magnetized neutron star whose cy- 
clotron line energy is positively correlated with luminosity 
dStaubertet al.l I2007t iKlochkov etafl 120121) . A counter ex- 
ample is V 0332+53, in which the cyclotron line energy in- 
creases as the source's luminosity decays and vice versa 
(Mowl avi et al.ll2006t iTsygankov et al.fe OloH. It was sug- 
gested by Staub ert et al.l d2007l) that apositive correlation be- 
tween E cyc and luminosity would occur for systems accreting 
at a sub-Eddington rate, while the negative correlation would 
occur when the accretion is locally super-Eddington. They 
show that the fractional change in E cyc is directly proportional 
to the fractional change in L. If such a relationship could be 
confirmed with cyclotron line sources of known distances it 
would be possible to use observations of cyclotron lines as 
standard candles. 

The 90% errors on the cyclotron line energy in Table [2] 
point to a possible correlation between £ cyc and luminosity. 

15 4U01 15+63 i s currently also conside red to be an example for a negative 
-E cyc -Z. correlation (Nakaiima et al. 2006). In this case the presence of the 
correlation, however, depends on the continuum model i Muller et al. 201 1). 



We calculated the 3cr errors on the line energy of LOl, L03, 
and DS3, observations which are clearly separated in time, 
and found that the energies are not consistent with each other 
within the errors. The 3cr confidence range on the CRSF en- 
ergy is 10.15-10.55keV during LOl, 9.58-9.85keV during 
L03, and 9.32-9.56 keV during DS3. The measured decrease 
in the line energy with decreasing flux (luminosity) is indica- 
tive of a positive correlation between E cyc and L and is appar- 
ent at the 99.7% confidence level in our data. Further obser- 
vations during a future outburst are necessary to confirm this 

c orrelation. 

Becker et al. (2012) recently derived a new expression for 
the critical luminosity assuming a simple physical model for 
the accretion column. Accordi ng to their estimate, and us- 
ing the B field determined in $4.11 the critical luminosity 
for Swift J1626.5-5156 is L crit = 1.5 x lO^flJ^ergs -1 
= 1.7 x 10 37 ergs~'. The normalized 3-20 keV luminosity 
for LOl, the observation with the highest flux F (Table |2), 
is L/L C rit = And 2 F IL aA = 1.6 at a distance d = 10.7 kpc, with 
all observations from DS1 and later falling below L cl i t . Thus 
Swift J1626. 6-5156 was mainly observed in the sub-critical 
accretion regime and the observed positive E cyc -L correlation 
is consistent with the theoretical expectation. 

4.4. CRSF parameter dependence on pulse phase 

The strong >20% variation of the fundamental energy 
with pulse phase that we measured for Swift J1626. 6-5156 
is c onsistent wit h sever al other sources, na mely Cen X- 
3 dBurderi et al.) 120001: ISuchv et all 120081). V ela X- l 
(lLa Barbera et all 120031 IKrevkenbohmetaLl 120021) . 
4U 01 15+63 (IHeindl et all 120041). and G X 301-2 
dKrevkenbohm et all 120041: ISuchv et all HoH. In all 
cases the fundamental energy varies by about 10-30% 
over the pulse. This variability has been attributed to the 
change in viewing angle and accretion column throughout 
the pulse, resulting in a variable local magnetic field as a 
function of pulse phase. Specifically the E cyc versus phase 
profile of Cen X-3 has very similar properties to the one of 
Swift J1626. 6-5156 (though shifted in phase from the falling 
to the rising edge of the pulse): a sharp rise followed by a 
slower decay. For Cen X-3 these E oyc variations were shown 
to be inconsistent with a changing viewing angle on a pure 
dipole magnetic field but it was found that it is possible that 
emission from above both poles is observed (Suchvetal. 
2008). The width and depth of the line are also generally 
variable over pulse phase, similar to what was observed for 
Swift J1626. 6-5156. We interpret the phase dependence of 
the CRSF parameters as additional evidence that the feature 
is real and not an artifact of the spectral fitting. 

4.5. Line spacing 

The ratio between the fundamental and harmonic line en- 
ergies for the phase averaged fit is 1.9 and ranges from 1.7 
to 2.0 for the four phase bins where the harmonic line is de- 
tected. Thus the harmonic probably does not occur at strictly 
twice the funda mental energy. Other sources, for exam- 
ple V 0332+53 dPottschmidt et al.1 120051: iKrevkenbohm et al 
120051: INakaiima et al.1 120101) and 4U 01 15+63 (IHeindl et al. 
ll999bt lSantangel o et all 119991: INakaiima et"afl 120061) have 
also displayed cyclotron lines at non-integer multiples of the 
fundamental energy. Spacing that is narrower than in the har- 
monic case is expe cted when gener al relativistic effects are 
taken into account (Meszaros 1992), while a wider spacing 
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has, e.g., been explained as being due the emission regions of 
the two lines being located at diff erent heights in the accre- 
tion column (Nakaii rriretaill2010l) . For Swift J1626. 6-5156 
the relativistic effect might have bee n observed, similar 
to V 0332+53 (Pottschmidtetal. 2005; Krevkenbohm et al. 
[2005 PI 

4.6. Iron line 

Swift J1626. 6-5156 is located in the direction of the Galac- 
tic plane, so the diffuse Galactic Fe K emission must at least 
contribute to our measured iron line flux. From the quies- 
cence o bserva tions we already know that this contribution is 
small ( !33.2.31 >. In addition we can make the following esti- 
mate: For a typical Galactic ridge location, the total Fe K 
line fl ux is ~ 3.8 x 10~ 4 photons s -1 crrT 2 deg~ 2 (Ebisawa et al. 
120081) . The PCA solid opening angle is 0.975 deg 2 and the 
effective area of one PCU at 6.4 keV is 1000 cm 2 , so the dif- 
fuse Fe K line flux (a veraged over the ar ea) seen by the PCA 
is ~0. 37 photons s -1 ( Jah oda et al . 1 120061) . We calculated the 
3-20 keV Fe K line flux from each of the three long observa- 
tions by first finding the flux of the best-fit model and then 
removing the emission line and, without refitting, calculating 
the flux again. We find Fe K fluxes of ~3.0 photons s _1 from 
LOl, ~2. 7 photons s _1 from L02, and ~1.0photonss _1 from 
L03, indicating that most of the iron line emission is coming 
from Swift J1626. 6-5156. We further note that the flux of the 
iron line, Ap e , as well as of the absorption A^h decrease over 
the time of the outburst (Table [2j. Since L oc M and since 
M depends on the amount of material in the system, this may 
reflect the diminishment of material from the neutron star's 
surroundings. 

5. CONCLUSION 

We have clearly detected a spectral feature at ~10keV in 
the spectrum of the Be/X-ray binary Swift J1626. 6-5156. Ev- 
ery aspect of our analysis points to this feature being a cy- 
clotron resonance scattering feature. The parameters of the cy- 
clotron line vary with pulse phase in a manner consistent with 
other accreting pulsars displaying cyclotron lines. The sec- 
ond harmonic of the CRSF, at an energy of ~18keV, is seen 
between pulse phases 0.4-0.7. The fundamental cyclotron 
line and the ~6.4keV iron emission line persist throughout 
the source's outburst and oscillatory stages. 

The neutron star's magnetic field strength of ~10 12 G, de- 
rived from the cyclotron line energy, is physically realistic 
and in the expected range for high-mass X-ray binaries. In- 
terestingly it is among the lowest B fields measured for an 
accreting pulsar so fa r, together with that of 4U 01 15+63 
(iNakajima et all 120061) . We find that this field strength is 
roughly consistent with the scenario of spin equilibration via 
disk accretion and we do not rule out the possibility that 
Swift J1626. 6-5156 may host an accretion disk. 

Comparing the 3<x confidence intervals of the fundamental 
CRSF energies in LOl, L03, and DS3 reveals a positive cor- 
relation between the line energy and the source luminosity. 
This behavior is similar to that of Her X-l and GX 304-1, in 
which the luminosity and line energy a re also positively cor - 
rela ted, and acco r ding to the model of Sta ubert et alj (2007) 
and iBecker etail (12012b implies that the outburst decay of 

16 For V 0332+53 the situation i s not c lear, however, since ratios >2 have 
also been obtained (Nakaiima et al. 2010), depending on the line and contin- 
uum model . 



Swift J1626. 6-5156 was sub-Eddington. Observations dur- 
ing future outbursts of Swift J1626. 6-5156 will be beneficial 
in confirming this relationship between its luminosity and cy- 
clotron line energy. 
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